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Abstract 

In this paper we derive analytical relations between probabilities of the excited state transfers 
and entanglements calculated by both the Wootters and positive partial transpose (PPT) criteria 
for the arbitrary spin system with single excited spin in the external magnetic field and Hamil- 
tonian commuting with I z . We apply these relations to study the arbitrary state transfers and 
entanglements in the simple systems of nuclear spins having two- and three-dimensional geomet- 
rical configurations with XXZ Hamiltonian. It is shown that High-Probability State Transfers 
(HPSTs) are possible among all four nodes placed in the corners of the rectangle with the proper 
ratio of sides as well as among all eight nodes placed in the corners of the parallelepiped with the 
proper ratio of sides. Entanglements responsible for these HPSTs have been identified. 
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I. INTRODUCTION 

This paper is devoted to the problem of the high probability state transfer (HPST) [l| 
among many nodes of the spin system and to the relationship between probabilities of HPSTs 
and entanglements responsible for these transfers. We consider nuclear spin-1/2 systems with 
XXZ-Hamiltonian and different geometrical configurations. However, the above relationship 
remains valid for any Hamiltonian 7i commuting with the total projection operator I z . 

WW 

By "state transfer" we mean the following phenomenon [2J, yj. Consider the chain of spin- 
1/2 with dipole-dipole interaction in the strong external magnetic field. Let all spins be di- 
rected along the external magnetic field except the ith one whose initial state is arbitrary. In 
other words, let the spin system be prepared in the state ipa = cos(0/2)|O) +e~ l ^sm(0/2)|z), 
where |0) is a ground state, i.e. all spins are directed along the magnetic field and |i) 
means that only ith spin is directed opposite to the external magnetic field (i.e. ith 
spin is excited). Let the energy of the ground state be zero. If the state becomes 
ipij = cos(6 l /2)|0) + e~ l ^fij sin(#/2)|j) with = 1 at the time moment t = t i3 - then 
we say that the initial state has been transferred from the ith to the jth node with the 
phase shift F^- = arg/y. Since \fij\ = 1, all other spins are directed along the field at t — tij. 
Here /y = fij(tij) and fij(t) is the transition amplitude of an excited state from the ith to 
the jih node: fij(t) = (j\e~ int \i) . It is known |2j that the effectiveness of the state transfer 
between the ith and jih nodes may be characterised by the fidelity Fy (t) 

= |^(t)|cosr tJ (t) |^(t)| 2 1 

tj U 3 6 2' 1 ' 

We see that the fidelity is maximal for = mod 2n. If the external magnetic field 
is homogeneous and we are interested in the state propagation between two nodes, say 
between sth and rth nodes at the moment t rs , then condition T rs = T rs (trs) = may be 
simply satisfied by the proper choice of the constant magnetic field value 2j. In this case 
the fidelity F rs = F rs (t rs ) takes maximal value together with absolute value of the transition 
amplitude |/ rs | = \f rs (t rs )\. For this reason, namely |/ rs | (rather then ^W)is considered as 
the characteristic of the state transfer in many refs, see, for instance, 0, 4, 5] 
| may 
P = I f I 2 

1 rs \JTs\ 

More general case of the HPSTs among many nodes of the iV node spin chain has been 
studied in jf|. In this case any particular state transfer between the ith and jth nodes is 



5| . It is clear that 

3e replaced by the probability of the excited state transfer P rs (t) = |/ rs (t)| 2 and 

13. 



associated with its own phase shift Tj,-, i, j 



, iV (note that i may be equal to j which 



means return of the state to the ith node). However, it is important that all these shifts 
may be eliminated using magnetic field properly depending on time [6J . For this reason the 
effectiveness of the state transfer between the zth and the jth nodes may be equivalently 
described either by the fidelity or by the probability of the excited state transfer Py even 
in this generalized case. Namely optimization of Py- allows us to find all necessary parameters 
of the geometrical spin configuration providing HPSTs among many nodes while phases 
may be removed by the appropriate time dependent external magnetic field as it was done in 
6|. For this reason, we will study the probability of the single excited state transfer instead 
of the fidelity (PQ) of arbitrary state transfer in the subsequent sections of this paper. This 
means that the spin 1/2 system of N nodes is prepared in the initial state 



(2) 



where i takes one of the values i — 1, . . . , N. 

The problem of the perfect state transfers (PSTs), HPSTs and entanglements in the spin 
systems is very attractive and different aspects of this problem have been studied in many 



details 



1 3 h, s s q, 



ii 
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14. Ilol. Nevertheless, most of the results are devoted 



to the linear and circular spin chains, which are considered as communication channels in 
the quantum information systems. Different Hamiltonians describing these chains have been 
studied, such as XY, XYZ, Heisenderg Hamiltonians. Usually, the approximation of the 
nearest neighbour interactions is taken as a basic tool for such studies. Note, that this 
is a good approximation in the case of, for instance, exchange interaction, when coupling 
constants decrease exponentially with increase in the distance. However, this approximation 
is not satisfactory for the spin systems with dipole-dipole interaction (such as nuclear spin 
systems in solids) and a wide spread of the coupling constants. 

Most efforts have been turned to the study of two phenomena: the state transfer along 
chains , llll, ll3] and two-qubit entanglements in chains (such as entangle- 



ment between end nodes or between nearest nodes {2, 



Il5l|). The Wootters criterion 



is applicable in this case allowing one to describe the entanglement in terms of so-called 



concurrence 



22j. It is important that there is an analytical dependence of the concurrence 
between two nodes on the probability of the state transfer between these nodes which was 

nn 

derived in |4J, [17j for the system with single excited spin. However, complicated system of 



N spins exhibits entanglements not only between two nodes, but also between arbitrary 



two subsystems. These entane 
transpose (PPT) criterion [18l . 



ements may be effectively described by the positive partial 
I9I introducing so-called double negativity as a measure of 
entanglement. Explicit relations of the double negativity associated with two subsystems 
on the probabilities of the excited state transfers between different nodes will be derived in 
this paper for the spin system with single excitation. 

Although the PST would be preferable in the quantum communication chains, it is hardly 
realizable in experiments with long chains because of the following two basic reasons: 

1. Theoretical prediction of the PST in the long chains is associated with the approxi- 
mation of the Hamiltonian by the nearest neighbour interaction, while the complete 
Hamiltonian must be used in practice. As we have already noted, this approximation 
is well applicable to the systems with exchange interaction and is not always valid for 
the systems with dipole-dipole interaction. 

2. Coupling constants may not be always known as accurately as we want in the case of 
both exchange and dipole-dipole interaction. 

Thus, HPST between different nodes of the spin system [lj seems to be more realistic in 
comparison with PST. 

In this paper we will study the spin systems with dipole-dipole interaction in the external 
magnetic field with single excited node described by the XX Z Hamiltonian. We will study 
the spin systems with different geometrical configurations of nodes and initial state (j2J) which 
may provide the HPSTs of the excited state among all of nodes. We will show which parts of 
the spin system must be entangled in order to provide each of these transfers. The study of 
two- and three-dimensional spin systems is important because they are more compact and 
consequently they are more promising as quantum registers and/or short communication 
channels. It will be shown in this paper, that namely such configurations (more precisely, 
spin configurations with nodes placed in the corners of either rectangle or parallelepiped) 
provide the HPSTs among several different nodes of the spin system during relatively short 
time interval in comparison with the line systems [6] which is important for the development 
of the quantum information systems and/or short communication channels. Our study is 



also stimulated by the experiments on the quantum information processes in solids 



24 



25] 
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This paper is organized as follows. In SecJTTl (and in Appendices [A] and [B]) we obtain 
analytical dependence of the either concurrences (Wootters criterion) or double negativities 
(PPT criterion) between different two subsystems of the spi n system on the probabilities 
of the state transfers, generalizing the results of refs. % ll7j] - In Sec JIVI we consider the 



simplest one-dimensional model of two nodes where the relationship between entanglement 
and probability of the state transfer is most transparent and an equivalent result may be 
obtained using either the Wootters |4j or PPT criterion. Two-dimensional spin systems will 
be considered in SecJV] see also Appendix O We arrange HPST among all nodes of the 
four-node spin system (rectangular geometry) and show that the external magnetic field 
directed along one of the sides of the rectangle decreases significantly (more then twice) the 
time intervals needed for the HPSTs among nodes in comparison with the case when the 
field is perpendicular to the plane of the rectangle. Similar study of the three-dimensional 
eight-node system (with spins placed in the corners of the parallelepiped) is represented in 
Sec JVIt see also Appendix [D] It is evident that HPSTs may not be effectively arranged in 
the arbitrary system of nodes. Detailed algorithm allowing one to obtain parameters of the 
rectangle spin system (namely, the ratio of sides of the rectangle) with the HPSTs among all 
four nodes is given in Appendix O Particular example of the eight-node three-dimensional 
spin system with HPSTs among all nodes (parallelepiped configuration) is represented in 
Appendix [Dl 

II. SPIN-1/2 SYSTEMS WITH SINGLE EXCITED NODE AND XX Z HAMILTO- 
NIAN 

We study the HPSTs and entanglements among the nodes of the spin-1/2 system in the 
external magnetic field described by the XX Z Hamiltonian with zero Larmor frequencies: 

TV 

^/ Dij(Ii tX l j t x Ii,ylj,y 2 1%.z I j,z) : (3) 

Dy = i - 3< r 52 S 2 «, (4) 

where 7 is gyromagnetic ratio, is the distance between zth and jth spins, 0y is the angle 
between the external magnetic field and r^, is the projection operator of the zth spin 
on the a axis, a = x,y, z, are the dipole-dipole coupling constants. This Hamiltonian 



describes the secular part of the dipole-dipole interaction in the strong external magnetic 
field 0- 

We denote D n = D nn+ i, n = 1, . . . , N — 1. Taking into account the definition of 
.Djj, for description of the spin system with arbitrary geometrical configuration we use the 
coordinates of each node, i.e. the set of the following triads 

(xi,yi,Zi), i = l,...,N, (5) 

so that 

Tij = yj {xj - Xi) 2 + (y,j - yi) 2 + (zj - z { ) 2 , (6) 

2 n ( Z 3 ~ z i) 2 

COS 9ij = — — 2 — — • 

It is important, that the Hamiltonian ([3]) commutes with I z (z-projection of the total spin): 

[H,I z ]=0. (7) 

This means that both 7i and I z have the common set of eigenvectors. It is convenient to 
write the eigenvectors of the operator I z in terms of the Dirac notations. Let 

|m...ntf) (8) 

be the eigenvector of the operator I z where the zth spin is directed opposite to the external 
magnetic field if Hi — 1 and along the field if = 0. For the sake of brevity, hereafter we will 
use notations |0) for the eigenvector associated with the state when all spins are directed 
along the external field and \ii . . .ik) for the eigenvectors associated with the state when 
iith, . . ., Zfeth spins are directed opposite to the external field, i.e. these spins are excited. 
Thus eigenvector \i) means that only iih spin is excited. Using these notations, the basis 
dHJ) may be ordered as follows: 

|»i>, |0), IM2), k < *2 < ■ ■ ■ < *jv, i-e. (9) 

ii = l,...,iV, i fc = ifc-i + 1, . . . , iV, k = 2,...,N. 

The matrix representation H of the Hamiltonian TC in basis ([9]) gets the following diagonal 
block structure: 

H = diag(if 1 , H , H 2 , H 3 . . . , H N ), (10) 
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where the block Hi is associated with the set of states of the whole spin system having i 
spins directed opposite to the field. 

Hereafter we will study the problem of the single excited quantum state transfer among 
nodes of the spin-1/2 system with the XXZ Hamiltonian in the external magnetic field. We 
say, that the /coth node is excited initially. It is important that only the block Hi is nonzero 
in this case: 



Hl = l(D-ri), 



(11) 
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where I is N x N identity matrix. This simplification of the Hamiltonian allows one to 

1. derive explicit analytical dependence of concurrence and/or double negativity (as mea- 
sures of entanglement between any two subsystems of the spin system) on the prob- 
abilities of the state transfers between different nodes of the system, which is hardly 
realizable in the case of the Hamiltonian with general structure (jTUJ) ; 



2. perform the numerical simulations of the state transfers in the big spin systems, which 



is hardly realizable in general case 



2l| 



Hereafter we will use the dimensionless time r, coupling constants d nm and distances £ nm , 



r = D i2 t, d 



1 wnm 



12) 



D12 r l2 

Using definitions (T121) and taking into account that D12 = ^h/r^ in all our examples, 
Hamiltonian ([3]) may be written as follows: 



A' 



Ti. — D^Ti., 7~L — dij(Ii iX Ij ;X + h, y Ij,y — 2Ii )Z Ij >z ), 



(13) 



III. PROBABILITIES OF THE STATE TRANSFERS AND ENTANGLEMENTS 
BETWEEN DIFFERENT SUBSYSTEMS OF THE ARBITRARY NUCLEAR SPIN 
SYSTEM 



First of all, in order to establish the relationships between probabilities of the state trans- 
fers among nodes and entanglements among them we will need density matrix p introduced 
as follows: 



P 



-iHt\ 



k ){k \e 



iHt 



(15) 



where k means that fcoth spin was directed opposite to the field initially. Only such initial 
states will be considered hereafter. It is important, that this matrix may be written in the 
following block form using basis 



( 



P 



A ®N,2 N -N 
®2 N -N,N ®2 N -N,2 N -N 



A 



a lN 



am ■ ■ ■ (Inn 



\ 



(16) 



where 0„ m means n x m zero matrix and nonzero elements defined as follows: 



,N, 



(17) 



a ij — (^|p|j) — fk ifk ji a ji ~ a iji h3 

where f nm are transmission amplitudes, 

N 

fnm (771 1 e \n) ^ 1 U n jU m jC J ' 

3=1 

Here Uij, 1, . . . , N, are the components of the normalized eigenvector Uj corresponding 

to the eigenvalue Xj of the matrix D: Duj = XjUj. 



fnm fmn • 



(18) 



A. Probability of the state transfer from the fcoth to the mth node of the A^-node 
spin chain 



The probability Pk om of the state transfer from the fcoth to the mth node as a function 

of time is defined by the diagonal element a mm of the density matrix. In fact 11], 

2 



Pk m(r) = (m\e- inT \k ) = \f 



Tkgm I 



(19) 
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Throughout this paper we will use notations Pk om and 7fc om for the probability of the 
HPST and for the time interval required for the HPST from the fcoth to the mth node of 
the iV-node chain: 

Pcom = Pkomij'kom) • (20) 

By our definition, the state transfer from the fcoth node to the mth node will be referred to 
as HPST if 

Pk m > Po- (21) 

The value Po is conventional, in our paper we take Po = 0.9. In addition, there is an 
important parameter of the HPSTs, namely the time interval T during which the excited 
state may be detected in all nodes of the system 

T = max r kon . (22) 

n=l,...,N 

B. Wootters criterion: two-node entanglements in the spin system 

It is well known that the entanglement between two nodes % and j of the iV-node spin 



system may be described by the Wootters criterion 
concurrence as a measure of the entanglement: 



221 ]. which introduces the so-called 



Cjj — max | 0, 2A — A n J , A = max(Ai, A 2 , A 3 , A 4 ). (23) 



ra=l 



Here A n , n = 1, 2, 3, 4 are the square roots of the eigenvalues of the 4x4 matrix 

Pij = (24) 
where * means complex conjugate, p\^ d is the reduced density matrix, i.e. 



p\? = \v, , ,,>. (25) 



Matrix p\ ed is defined as 



Pt = °1® °W° V i ® (26) 



After simple calculations (see Appendix |A] for details) one derives the following formula: 



This relation is valid 
commuting with I z 0, ~v\ . 



Cij — 2\a,ij\ — 2a/ PcoiPfcoj) i 7^ 3- (27) 
or the system with any number of spins and for any Hamiltonian 
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C. PPT criterion: entanglement between two arbitrary subsystems 



PPT criterion describes the entanglement between any two subsystems A and B of the 
system S in terms of the so-called double negativity Na,b jlOl, |23j, which is the absolute 
value of the doubled sum of the negative eigenvalues of the matrix pab-,c, 

Pab-,c = (pa,b) Ta , S = AUBUC, (28) 

where the reduced density matrix pa,b is defined as follows: 

Pa,b = Tr c p, (29) 

and Ta means the transposition with respect to the subsystem A. In particular, the subsys- 
tem C may be empty. It is important that one can write the explicit formulae for Na,b f°r 
the spin system in the magnetic field with single excited spin described by any Hamiltonian 
commuting with the total projection I z (see Appendix [B] for details): 



N, 



,3i---3m 2 



<?h...i 



m 1 3i---3m 2 



Ah M 2 



(30) 



n=l m=l 



a 



H.-%M 1 3l---3M 2 



Ah 



M 2 



■■ 1 m 1 3i---3m 2 



where A = {i 1 . . . i Ml }, B = {j x . . . j M2 }, 



a 



n...ijsr n 



we use the identity J2^=i P fcon = !)• ^ n particular, 



N 

E 



N 



n=l 



No 



m=l 



(31) 



n=l 



n=l 



n=l 



h3 



i.l 



Or 



a ij + 4Pfc jPfc j 



(32) 



^1 - P ko i - Pk j - \J0- - Pk i - Pk j) 2 + 4P fcoi Pfc j^ 



i,rest 



N 



\ 



N 



Pk i z2 Pk °i ~ 2a / Pkoii^ - Pkoi)-, 



3 = 1 



N, 



W2,3l32 



Hi23l32 



\ 



2 2 



(33) 



(34) 



n=l m=l 



a 



^hnhh ^{Pkoh + Pkni2){Pkoji + ^Wta)) 
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where Ni^ rest is double negativity associated with the entanglement between ith spin and 
the rest part of the spin system. It has the most simple expression depending only on Pk i- 
Eqs. (I30|321l34l) are simple relations between two important characteristics of the spin 
system: entanglement between arbitrary subsystems A and B (characterized by the double 
negativity) and probabilities of the state transfers from the fcoth node to other nodes of the 
subsystems A and B. However, since (a) both these characteristics have rather complicated 
oscillating behaviour and (b) it is difficult to "inverse" relation ( l30i) (i.e. to write probabilities 
as functions of double negativities), it is hard to get answers to the following questions: 

1. Which geometrical configurations of the spin system provide HPSTs among all nodes? 

2. Which entanglements are responsible for the HPSTs among different nodes of the spin 
system? 

We will answer these questions for the particular systems of two, four and eight nodes in 
SecsfjVHVI] 

IV. PERFECT STATE TRANSFER AND ENTANGLEMENT IN THE SYSTEM 
OF TWO NODES 

n 

It was shown in [3J] that the probabilities for the excited state to be detected in the first 
and in the second nodes are following (ko = 1): 

P n (r) = cos 2 (r/2), P 12 (r) = sin 2 (r/2), (35) 

while the double negativity iVi )2 is defined by eq.(l32l): 

Ni, 2 (t) = C 12 (r) = 2y/Pn(r)P 12 (r) = | sinr|, (36) 

i.e. P\2 = 1 at Ti2 = 7r + 2-nn, n = 0, ±1, ±2, . . .. We compare the graphs of -Pii(t) and 
Pi 2 (r) with double negativity AT 1)2 (r), see FigHJ and observe that maxima of iVi >2 correspond 
to Pu = P\2 = 1/2 while the minima of iV lj2 correspond to the maxima and minima of Pu, 
% = 1,2. This result is not surprising. In fact, the wave function associated with the two- node 
spin system is following: 

*(r) = /u(r)|10) + / la (r)|01) J P n = \f n \ 2 , Pi2 = \fn\ 2 . (37) 



11 




FIG. 1: The probabilities P\\ and P\2 and the double negativity N\ 2 in the two-node spin system 

However, if P u = or P 12 = this function reduces to \& = 1 01) or ^ = 1 10) respectively. It 
is well known that these states are separable and their entanglements are zero. Vice-versa, 
if Pu = P\2 = 1/2, then we have the sing let state * = l/\/2(|01) - |10)) which is the most 
entangled one. 

We see from FigJT] that the relation between the probability of the state transfer and 
the entanglement in the two-node system is very transparent. Moreover, we will see that 
relations between probabilities of state transfers and entanglements in more complicated 
systems are very similar. However, the system with N > 2 requires additional analysis in 
order to find such geometrical configuration of the spin system which allows the HPSTs 
among all nodes (or, may be, among some of them). Algorithms allowing us to perform this 
analysis for the rectangular system of four nodes and for the eight-node system with spins 
placed in the corners of the parallelepiped are represented in Appendices O and [D] 



V. HPSTS AND ENTANGLEMENTS IN THE SYSTEM OF FOUR NODES WITH 
RECTANGULAR GEOMETRY 

We consider the rectangle shown in Figj2j Here £12 = 1 while the length of another side 
is a parameter of the system, £ 14 = b. However, it is more convenient for us to introduce 
parameter 



12 



2 



1 *- 



ft 4 =b 



FIG. 2: The rectangular system of four nodes. Here £12 = 1 while £14 is a parameter: £14 = b. We 
will use parameter 5 = 1/6 3 instead of b. 

instead of b. Due to the symmetry of the spin system, matrix D in eq. fllll) ) reads: 

/ 2f 1 d 13 d u \ 
1 2T di4 di3 
d 13 du 2T 1 
\ d u d 13 1 2f / 

The structure of the matrix D helps us to guess the following set of independent normalized 
eigenvectors: 



D 



l + d 13 + d 



114. 



(39) 



u x = 1/2(1 -1 1 -If, u 2 = 1/2(1 1 1 If, 

m 3 = 1/2(1 - 1 - 1 if, m 4 = 1/2(1 1 - 1 - if 



(40) 



with the appropriate eigenvalues: 



Ai = 2T - 1 - d 1A + dia, A 2 = 2r + 1 + d u + d 



13; 



A 3 = 2r - 1 + d 14 - d 13 , 



2T + 1 - d u - d 13 . 



(41) 



It is remarkable that 

1. these eigenvectors do not depend on b and 

2. any component of any eigenvector is either 1 or -1. 

Namely the latter property guaranties the HPSTs among all nodes of the four-node system. 

Due to the symmetry of the rectangular cluster, it is enough to consider the case with 
initial excited state in the first node, i.e. k = 1. Eq.( jT9i) yields explicit expressions for the 
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probabilities P\f 

Pn = -^1 + cosr^cos(<ii4r) + cos(<ii3r)j + cos(<ii4r) cos(<ii3r)j, (42) 

P12 = - ^1 — cosr^ cos(c?i4r) + cos(c?i 3 r) j + cos(g?i4t) cos(d 1 3r) j , 

Pi3 = -^1 + cosr^cos(<ii4r) — cos(<ii3r)j — cos(rfi 4 r) cos(rfi3r)j, 

P$ = - (l + cos t ^ — cos(di4r) + cos(c?i 3 r)^ — cos(di4r) cos(di 3 r) 



The explicit formulae relating the coupling constants dij and b depend on the direction of 
the external magnetic field in accordance with the definition of , see eq. (1141) . We consider 
two following cases: 

1. The external magnetic field is perpendicular to the rectangle, so that 

d u = 5, rf 13 = (l + ft 2) 3/2 - ( 43 ) 

2. The external magnetic field is directed along the side b of the rectangle, so that 

d u = -25, ^3 = (l " yf^) (1 + &T 3/2 - (44) 

Remark that both cases ( ]4"3l) and ( 1441 allow equality {d^l = 1 if b = bo = 1 for case ()4*3l) 
and b = b = 2 1 / 3 for case ( ]44"1) . Rectangles with these special values of b may not provide 
HPSTs to all nodes. In fact, the appropriate expressions for Py become simpler in this case: 

P n = - (l + cos 2 r + 2cosrcos(rfi 3 r)) , (45) 

Pn = Pu = ^-p < \ < Po = 0.9, 
Pi3 = - (l + cos 2 r — 2 cost cos (g? 13 t)) . 

We see that in this case maxPi2 = maxPi4 = 1/4, i.e. HPST exists only between 1st and 
3rd nodes. Any rectangular configuration with b ^ bo provides HPSTs among all nodes, 
although the appropriate time interval T may be long. 



A. Relationships between probabilities and double negativities 



As we noted above, any rectangle with b ^ bo provides HPST among all nodes. However, 
in general, the time interval T is long. An important problem is to find such values of b 

14 



which provide relatively short time interval T. This problem is solved in Appendix O for 
the rectangular configuration of four nodes, where we have found the values of b such that 
T < 10 and T < 15 for the case with magnetic field perpendicular to the plane of rectangle: 

5 G [5.56, 9.62] for T < 10, (46) 
5 E [5.56, 17.79] for T < 15. 

and values of b such that T < 3.5 and T < 6 for the case with magnetic field directed along 
the side b: 

5 E [2.62, 6.08] for T < 3.5, (47) 
5 E [2.32, 6.08] U [14.89, 30.63] for T < 6. 

Thus we see that T is longer in the first case (magnetic field perpendicular to the plane of 
the rectangle), i.e. the second case is more preferable for the organization of the HPSTs. 

To demonstrate the qualitative relationship between the probabilities of the state transfers 
and the double negativities (see eqs. (l30f321l3"4"|) ) we show their graphs corresponding to the 
case (1441) and 5 = 4.3 in Figj3] during the interval T = [0, 3.5]. We see that the whole interval 
may be separated into three parts. During the first part, < r < 1.1, the probabilities P\\ 
and P14 have big amplitudes. The associated big amplitude double negativity is showing 
that the first and the 4th nodes are most entangled during this interval. During the last 
part, 2.2 < r < 3.5, the probabilities Pyi and P13 have big amplitudes. The appropriate big 
amplitude double negativity is jV 2) 3 showing that the second and the third nodes are most 
entangled during this interval. The middle part, 1.1 < r < 2.2, is characterized by the small 
values of P\j and considerable double negativity iV 14j 23 showing us that two opposite sides 
(namely sides 1-4 and 2-3) of the rectangle are well entangled. Two more double negativities, 
AT 1)2 and iV^, are also considerable during the middle interval but they are not represented 
in the figure because their role is equivalent to the role of ^14^3. The double negativities 
AT 1)3 and A^ 24 remain small during the whole interval T, see Figj3]Jc). We conclude that the 
time interval T = t 13 3.29. 
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T 



FIG. 3: Four- node system with the external field directed along £x4- (°0 The probabilities Pu, 
(ti4,Pi 4 ) = (0.36,0.97), (r^A) = (2.92,0.96), (m,P 13 ) = (3.29,0.96); (6) the double nega- 
tivities which provide HPSTs in the system; (c) the two-node double negativities which are not 
associated with HPSTs; 

VI. THREE-DIMENSIONAL SPIN SYSTEM OF EIGHT NODES WITH HPSTS 
AMONG ALL OF THEM 



In this section we consider the parallelepiped with spin-1/2 nodes in its corners, see FigJH 
Matrix D in eq. (fTT|) has the following block structure: 

^ R\ R2 R3 R4 ^ 
R2 R± R4 R3 

R3 R4 Ri R2 

\ Ra R3 R2 R\ ) 

2f 1 \ I d 13 d u 

1 2f J ' ' ' \ d 14 d 13 



D 



(48) 



Ri 



R* 




Ra 
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FIG. 4: The three-dimensional eight-node system. The external magnetic field h is directed along 
the side £15. Here £12 = 1, while two other sides are parameters of this configuration: £14 = b%, 
£15 = ^2- We will use parameters <5j = l/bf, i = 1, 2 instead of 61 and 62- 



du = (l + ^f 3/2 , d u = bf, d lb = -2b 



-3 
2 5 



^=( 1 - 3 ^T6f)^ + f '?)" 3/2 ' f = |>j. ^ = 1. 

Again, the structure of the matrix D allows us to find the following set of normalized 
independent eigenvectors: 

«i = mm if, «a = ^=(1 1 1 1 " 1 " 1 " 1 " if, (49) 

u 3 = ^=(1 1 - 1 - 1 1 1 - 1 - if, U4 = _L(1 1 - 1 - 1 - 1 - 1 1 if, 

u 5 = ^=(1 - 1 1 - 1 1 - 1 1 - If, Ue = ^=(1 - 1 1 - 1 - 1 1 - 1 If, 

u 7 = ^=(1 - 1 - 1 1 1 - 1 - 1 If, u 8 = ^=(1 - 1 - 1 1 - 1 1 1 - If. 

Similar to the case of four nodes, the eigenvectors do not depend on b i: i = 1,2. The 
correspondent set of eigenvalues reads: 

A 1= 3f, A 2 = 3f -2(d 15 + d w + d 17 + d 18 ), (50) 
A 3 = 3f - 2(d n + d u + d 17 + d 18 ), A 4 = 3f - 2(d 13 + d 14 + d 15 + d w ), 
A 5 = 3f - 2(1 + d u + die + d 18 ), A 6 = 3T - 2(1 + d u + d lb + d 17 ), 
A 7 = 3f - 2(1 + d 13 + d 16 + d 17 ), A 8 = 3f - 2(1 + d 13 + d 15 + d 18 ), 

Due to the symmetry of our cluster, it is enough to study the case with the initial excited 
state in the first node, i.e. k = 1 similar to SecJV] The expressions f[T^]) for the probabilities 
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Pij in terms of the coupling constants dij are rather complicated so that we do not represent 
them here. Note, however, that the cube does not allow the HPSTs among all nodes. In 
fact, in the case b± — hi — 1 one has 

Pl2 = Pl4 = ^ < 1 (51) 



16 ~ 16' 



Pyy = — ( 7 + cos(4r) + 8 cos — T -= ( cost + cos(2r) cos — T -= 1 + 
32 V 4^2 V 4^2, 

3r 

4 (cost + cos(3r)) cos 



4^2 



Pi 3 = — f 7 + cos(4r) — 8 cos — = (cost — cos(2r) cos — =| 
32 V K } 4^2\ K 1 AV2J 

3r 

4 (cost + cos(3t)) cos 



4V2 

sin 2 t f 3r \ 

Pis = I 3 + cos(2r) + 4 cost cos — -= I , 

8 V 4V2/ 

Pie = Pi8 = ^(3 + cos(4r)-2cos lV 4 ; -2cos lV ^ ' j < -, 
sin 2 r f „ ,„ x 3r \ 



Pi7 = ( 3 + cos(2t) — 4 cos t cos ■ _, 

8 V 4^7 

so that the probabilities Pi 2 , Pi 4 , Pi 6 and Pig may not approach P = 0.9. 

The study of the HPSTs and entanglements in this system requires finding such subspace 
of the two-dimensional space of positive parameters 8\ and 62 that any pair (5i,52) from 
this subspace provides the HPSTs among all nodes. This computational problem will not 
be considered in this paper in the full extend. Instead, we consider an example. Namely, 
we will show (see Appendix |D] for details) that the HPSTs among all nodes are possible for 
81 = 9 and 82 = 26.20 during the r-interval T < 25. 

We demonstrate that the relationship between probabilities and double negativities is 
very similar to one considered in Sec]V] For this purpose, let us refer to FigsJMH FigsJS^a) 
and (b) collect all probabilities and all double negativities involved into the state transfer 
process. Similar to the four-node system considered in Sec]V] we select three parts of the 
whole interval T = [0,251: < r < 0.5, 0.5 < r < 23 and 23 < r < 24. It is clear from 
FigE^a), that the HPSTs take place in the first and in the third parts of T. All probabilities 
of the state transfers are not high during the second part of the above interval. Amplitudes 
of the probabilities Ph, i = 1,4,5,8, are big during the first part of the interval, while 
amplitudes of the probabilities Ph, i = 2,3, 6, 7, are big during the third part of the interval 
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FIG. 5: Eight-node system. The probabilities and double negativities corresponding to S\ = 9, 
82 = S = 26.20. The HPSTs take place during the first and the last parts of the time interval 
< r < 25. 

T. Similarly, Figj5^&) shows that double negativities iVi s, N^ 8 and Ni 5 ^$ are significant 
during the first part of the interval T while double negativities N 2 fi, iV 3j7 and A^6,37 are 
significant during the third part of this interval. One more double negativity ^458^367 is 
significant during the second part of the interval T and is not high during the first and the 
last parts. This means that namely ^1453,2367 is responsible for the HPSTs from the plane 
1-4-5-8 to the plane 2-3-6-7. 

The probabilities and double negativities during the first and the third parts of the 
interval T are represented in FigsE] and respectively in more details. We show only those 
probabilities whose amplitudes exceed the value Pq = 0.9. One can see from FiglB^a, b) that 
the HPSTs to the 4th, 5th and 8th nodes occur during the interval < r < r 14 = 0.36. 
Functions and provide the HPSTs between the first and the 5th and between the 
4th and the 8th nodes respectively, while -/V 15; 4 8 provides the HPST from the side 1-5 to the 
side 4-8. 
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Similarly, we see from FigJ7|(a, b) that the HPSTs to the 2nd, 3rd, 6th and 7th nodes 
occur during the time interval (ti 6 = 23.23) < r < (ri 2 = 23.89). Functions N 2 ^ and jV 3) 7 
provide the HPSTs between the 2nd and the 6th and between the 3rd and the 7th nodes 
respectively, while N 26j37 provides HPST from the side 2-6 to the side 3-7. We also conclude 
that T = r 12 = 23.89. 




FIG. 6: Eight-node system. The probabilities and double negativities corresponding to 5\ = 9, 
5 2 = S = 26.20 and HPSTs during the time interval < r < 0.5; (t 15 ,P 15 ) = (0.06,0.93), 
im, Pi 8 ) = (0.30,0.91), (t u ,P u ) = (0.36,0.93) 




FIG. 7: Eight-node system. The probabilities and double negativities corresponding to 5± = 9, 
5 2 = 6 = 26.20 and HPSTs during the time interval 23.1 < r < 24; (t w ,P w ) = (23.23,0.90), 
(ri3,Pi 3 ) = (23.53,0.95), (t 17 ,P 17 ) = (23.59,0.95) , (r 12 ,P 12 ) = (23.89,0.91) 



20 



VII. CONCLUSIONS 



We have derived simple relations between the probabilities of the excited state transfers 
to different nodes of the spin system and the entanglements between different parts of this 
system described by PPT criterion for the arbitrary spin system in the external magnetic 
field with single excited spin and any Hamiltonian preserving the number of excitations, 
such as XXZ Hamiltonian. Although similar relations for the concurrence (which is a 
measure of the entanglement between two nodes in accordance to Wootters criterion) has 
been found 4], [l7|, the PPT criterion allows one to involve entanglements between two 
arbitrary subsystems of nodes which is important for the systems with HPSTs among many 
nodes as has been illustrated in this paper. 

We have found examples of 4- and 8-node spin systems which provide HPSTs among all 
nodes. We have seen that the HPSTs between two subsystems require the high entanglement 
between them. This is illustrated in all examples considered in this paper: two-node, four- 
node and eight-node nuclear spin systems, see Figsfj], [31 [5H3 respectively. 

In all spin systems considered in this paper, the HPSTs among many nodes are possible 
due to the remarkable property of the eigenvectors of the matrix D (see eq. ffTT]) ): all their 
elements are real, equal by absolute value and deffer only by sign. It is important that four- 
and eight-nodes systems with this property have simple geometrical configuration: either 
rectangle (four nodes) or parallelepiped (eight nodes). It is possible to construct the higher 
dimensional eigenvectors with this property using the following simple algorithm. 

Consider the eigenvector spaces with the above mentioned property: all elements of all 
eigenvectors are real and have the same absolute value. Suppose that we have M-dimensional 
space Bm of such vectors. Then we may construct 2M dimensional space by the formula 



& 



2M 



1/V2 ((1,1)®B m U(1,-1)®Bm). (52) 



Thus we are able to construct 2 s -dimensional basis with the above property. However it 
is difficult to find appropriate spin configurations different from those which have been 
described above. Of course, it is quite possible that the HPSTs among many nodes may 
be arranged using completely different mechanism allowing one, for instance, to handle the 
coupling constants in complicated spin systems 0, Q , which is one of the open problems 
for further study. 
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APPENDIX A: WOOTTERS CRITERION. DERIVATION OF EQ.fl27|) 



Our calculations are based on eqs. fl23H25|) . By construction, the reduced density matrix 



pred defined ^ eq. (!25p takes a simple 4 x 4-dimensional form in the following basis 

1 10), |01), |00), |11), 



(Al) 



where the first and the second elements are assotiated with the ith and the jth nodes 
respectively. One has 



( an ciij \ 



aji ajj 
Gij 

\ o o ooy 



N 

E 



N 



Pi 



hot 



Pi 



k j- 



n = l 



Then 



( ajj a^ 0^ 



o o 



(A2) 



(A3) 





V a iS J 

Direct calculation shows that the matrix (J24l) has only one non-zero eigenvalue: Ai 
4\aij\ 2 = 4P koi P koj , which yields eq.flSZJ). 



APPENDIX B: PPT CRITERION. DERIVATION OF EQ.([30]) 

We derive formulae fl30|) in this section. First of all, in order to calculate A^..,^ ji...^ 
one needs the reduced density matrix Pii...i M ,h-3M calculated in accordance with eq. (l29|) 
where 

A = {i u ...,i Ml }, B = {j u . . . , j M2 }. (Bl) 

This is 2 Ml+M2 x 2 RIl+M2 matrix. However, most of elements of this matrix are zeros. All 
nonzero elements are collected in the K x K (K = M\ + Mi + M\Mi + 1) block on the 
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diagonal of Pi 1 ...i Ml ,j 1 ...j M2 - This block corresponds to the subspace spanned by the following 
set of vectors (we use notations equivalent to ones introduced for basis 



\in), \jm), |0), \injm), U = 1, . . . , M X , 171 = 1, . . . , M 2 . 



(B2) 



We refer to this block as Pi 1 ...i Ml j 1 ...j M2 

( 



Ph. 



■lM 1 3l---3M 2 



Qii...iM 1 ji—3M 2 



Oi 



Hi 



(J h—iM x h—3M 2 



\ 



of, 



oh 







IV 



(B3) 



/ 



where Q^...^ j 1 ...j M is (Mi + M 2 ) x (Mi + M 2 ) square matrix 



Qh—iM 1 jl—3M 2 

( 



R T 



Ri 1 ...i Ml Rii-iM 1 h-3M 2 



(B4) 



n—iMi31---3M 2 



\ 



M X 



( 



Rjl-3 



3M 2 



V a «M 1 *l ' ' ' a iM 1 iM 1 J 

( „ 



"3i]i 



a 



313M 2 



Ri 1 ...i Ml j 1 . 



■3 M 9 



a il3M 2 



\ 



y a iM 1 ji ' ' ' a iM 1 3M 2 J 

3m 2 = ( a ikji cdots &i k j M ), 



1h;ji 



\ 9iM 1 ;jl-3M 2 f 



\ a 3M 2 3l ' ' ' a 3M 2 3M 2 J 

\ 



-3M 2 



0/ is the column of Mi + M 2 zeros, 0// is the (Mi + M 2 ) x MiM 2 zero matrix, 0/// is the 
row of MlM 2 zeros, Ojy is the MiM 2 x MiM 2 square matrix of zeros. 

Now we have to transpose the matrix Pi x ...i M j\-3m with respect to the nodes i ly . . . ,i Ml . 
The nonzero diagonal block of the resulting matrix is associated with transposition of the 
block Pii...i M ji...j M with respect to the nodes h, . . . ,iMi' 

l 31 ■ ■ -3 ai 2 * 



~"± •••"JWl 

Ph—iMi3l—3A 



\ 



0] 
0/j 



<7, 



l—i Ml h—3M 2 



R 







(B5) 



/ 



where + means hermitian conjugation, 



Qh...i Ml 3i- 



■3 m,, 



R ii-i Ml °v 

0y Rji~3M 2 



Ri — {Qii;ji—j 



3m 2 ■ ■ ■ (JWpjl— 3M 2 )i 



(B6) 
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Oy is the Mi x M 2 zero matrix. By construction, both blocks Rf^^ and Rj 1 ...j M2 are density 
matrices and, as a consequence, have non-negative eigenvalues. Thus, all eigenvalues of the 
matrix Q r - ed , „• „■ are non- negative. For this reason, looking for the negative eigenvalues 
of the matrix 5* 1 „• , we stay with the matrix 



<J h—iM 1 h-0M 2 ^ I 

Rf Oiv 



(B7) 



Its characteristics equation reads: 

(_ X) M lM2 -l | A 2 /W 



n=l m=l 



:bs) 



which has two nonzero roots one of which is negative: 



•^1 9 I a i\-iM 1 h- 



■3 Ah 



'l-'M, 31— 3 M* 



+ 



Mi M 2 
n=l m=l 



(B9) 



Consequently, 



2|Ai|, 



(BIO) 



which generates the formulae fl30l) . 



APPENDIX C: OPTIMIZATION OF THE RECTANGULAR SYSTEM. VALUES 
OF THE PARAMETER 5 PROVIDING THE HPSTS AMONG ALL NODES 

There are four functions characterizing HPSTs in the system of four spin- 1/2 nodes: 



PuM, 2 = 1,2,3,4. 



(CI) 



As for the entanglements, the situation is more complicated. We have the following lists of 
all possible double negativities in the four-node system: 

tf ltl = {N id , i,j = 1,2,3,4, i^j}, (C2) 
A/" 2 , 2 = {A^ 12 ,34, N 13i24 , N U:23 }, (C3) 

M,3 = {#1,234, #3,134, jV 3 ,124, #4,123}- (C4) 

We will show that Mi x and A/2,2 are responsible for the HPSTs among all nodes. 
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Our system has one geometrical parameter b or 5 (see eq. ( 1381) ) completely describing the 
rectangular geometry. We represent an algorithm allowing one to find such values of the 
parameter 5 which provide inequality (|2ip with the short time interval T fl22|) . For this 
purpose we consider the functions Pu and Nij as functions of two arguments, r and 8. 
Let us fix some time interval T and check whether the state may be transferred with high 
probability from the 1st node to any other node and whether the entanglements between 
any two nodes are significant during this time interval. For this purpose we construct two 
following functions: 

F P (5,T)= lim F p (5,T,Ar), (C5) 
F N {5,T) = \im F N (8,T, At), (C6) 



where 



fc=l,2,3,4 



i=0,...,K 



F p (8,T,At) = min max P lk (n, S) , (C7) 



F (S,T,At)= min max N nm (Ti,8) 

.i=0, ...,K 



n,m=l,2,3,4 



T 

n = iAr, At = — , KgN. 



The HPSTs among all nodes of the system are possible if there is such 8 = 8q that 

F P (8 ,T)>P . (C8) 

Then 

T <T. (C9) 

Function tells us how significant is the entanglement between any two nodes of the 
system. 

First, we consider the case when the external magnetic field is perpendicular to the 
rectangle. Functions F P (S,T, At) and F (8, T, At) with At = 0.01 corresponding to the 
intervals T = 10 and 15 are represented in FigJHJ The HPSTs among all nodes are possible 
for the parameter 5 inside of the intervals (|46|) . 

The interval of 5 corresponding to the HPSTs among all nodes increases with increase in 
T. If T is big enough (we have found that T > 40), then the HPSTs among all nodes exist 
even for S < 1. If 8 = 8q = 1, then F p may not exceed 1/4 in accordance with eqs. (1451) . 
FigE] shows that the function F N ~ 0.8 0.9 when F p > 0.9. This means, that the HPSTs 
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FIG. 8: Four-node system with the external field perpendicular to the plane of the rectangle. 
Functions F p (<5, T, 0.01) and F N (S, T, 0.01). (a) T = 10, $i = 5.56, 5 2 = 9.62; (6) T = 15, 
5 1 = 5.56, <5 2 = 17.79. 

among all nodes during the interval T entangle any two nodes in this case, i.e. functions 
Ni t j (for i,j = 1, 2, 3, 4, i ^ j) must have big amplitudes during the time interval T. 

Similarly, the functions F P (S, T, At) and F N (S, T, At) for the case with the magnetic field 
directed along £ i4 are represented in Figj9]for, T = 3.5, At = 0.01 and T = 6, At = 0.001. 
The HPSTs among all nodes are possible for 5 inside of the intervals (|47l) . In this case the 
function F N does not necessary take a big value when F p > 0.9. This means, that not any 
two nodes must be entangled in order to provide the HPSTs among all nodes during the 
time interval T. Note that exceptional value of S is S = 5q = 1/2, when F p does not exceed 
1/4 in accordance with eqs. (1451) . 




FIG. 9: Four-node system with the external field directed along £14. Functions F P (8,T) and 
F N (5,T); (a) T = 3.5, 6 X = 2.62, 5 2 = 6.08, S = 4.3. F p (5 0l 3.5) = 0.96 (6) T = 6, 6% = 2.32, 
6 2 = 6.08, 6 3 = 14.89, 64 = 30.63. 

Comparison of FigsJHta) and[9](a) shows that the second case (i.e. external field is directed 
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along £14) is more preferable for the organization of the HPSTs because appropriate interval 
T is almost three times shorter. 

FigsJHJ and [9] show that the intervals of the parameter 5 providing the HPSTs among all 
nodes increase with increase in T, i.e. the system becomes more "stable" with respect to 
variations in b (compare the interval 5i < 5 < 62 in Fig 18] and the intervals 5i < 5 < 82, 
S3 < 5 < 5 4 in FigE]). 

Finally, remark that function ( 1C7I) are decreasing functions of At. Thus, one may expect 
that the graphs of the functions F P (5,T) and F N (5,T) are above of the appropriate curves 
shown in FigsJH] and [9] However, further decrease in At negligibly effects the shapes of the 
above curves in all examples of this section. 



APPENDIX D: ON THE OPTIMIZATION OF THE EIGHT-NODE SYSTEM. 
VALUES OF THE PARAMETERS <5i AND 5 2 PROVIDING THE HPSTS AMONG 
ALL NODES 

We introduced functions F p and F N in Appendix [C] see eqs. (]C5IIC6l) . However, we have 
seen that F N is not as helpful as F p in defining the optimal parameter 8. In fact, FigsJHJ 
and [9] show us that F N may be either big or small when F p > 0.9 reflecting the fact that 
not all entanglements between two nodes are important for the HPSTs during the interval 
T. For this reason we consider only F p in this section: 



F P (5,T) = lim F p (8,T,At), 



;di) 



where 



F p (5, T, At) = min max Pik{T;,8) 

fc=l,..,8 li=0,...,K 



n = iAr, At = —, K e N, (D2) 
K 



and Pik are given by eqs. f)19ll49ll50|) . 

Function F P (5,T, At) is shown in Fig{]7j for T = 25 and At = 0.01. We see from 
Figfiorfc) that the HPSTs among all nodes are possible, for instance, for 5 = 5q = 26.20. 
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